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Abstract. The geometric theory of Lie systems wiU be used to estabhsh integrabihty 
conditions for several systems of differential equations, in particular, Riccati equations 
i-^ . and Ermakov systems. Many different integrability criteria in the literature will be 

analysed from this new perspective and some applications in physics will be given. 
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1. Introduction: Non-autonomous systems of differential equations 

'i> I Non-autonomous systems of differential equations play a relevant role in many physical 
^ I theories, describing evolution in terms of either time or any other parameter. For 
O ' instance Hamilton equations are systems of first-order differential equations and Euler- 
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Lagrange equations for regular Lagrangians are systems of second-order differential 
equations. These systems also appear in many other problems, and the theorem of 
local existence and uniqueness of solutions is of a crucial importance. However, it is not 
possible to finding exphcit solutions of such systems and approximation methods have 
^ been developed for dealing with them. This gives even more relevance to the explicitly 

^ ' solvable examples which can be used to test approximation methods. In the linear case 
there is a linear superposition principle allowing us to write the general solution as a 
linear combination of a fundamental set of solutions. Particular solutions are also used 
to reduce the problem to a simpler one (see e.g. [T] for Riccati equations). Our aim in 
this paper is to analyse sufficient conditions for integrability of a particular case of non 
autonomous systems, the so called Lie systems [21 El IH E] . We shall understand that to 
find a solution means to reduce the problem to some quadratures. 

We remark that a system of second-order differential equations in n variables is 
related with a system of first-order equations in 2n variables, 

x' = F\x,x,t) ^ \ ~ ^^ii , i = l,...n. (1) 

I ^ F\x,v,t) 

and therefore it is enough to restrict ourselves to study systems of first-order differential 
equations. For instance, for the harmonic oscillator with a t-dependent angular 
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frequency u(t): 



From the geometric viewpoint, a system 

x' = X'{x,t), i = l,...,n, (3) 
is associated with the t-dependent vector field 

X = X\x,t)^ (4) 
whose integral curves are determined by the solutions of ([3]). In particular, for ([2]), 

d 2/ \ 9 

X = V- u {t)x — , 

ox ov 

The general solution of the inhomogeneous linear differential equation 

'^ = ho{t) + hi{t)x 
can be found with two quadratures. It is given by 

= exp hi{s)ds^ x ^Xq + j bo (t') exp 

and when the systems we are dealing with are linear, there is a linear superposition 
principle allowing us to find the general solution as a linear combination of n particular 
solutions. For instance, for ([2]) if we know a particular solution the general solution 
can be found by means of one quadrature, while if we know two particular solutions, 
Xi and X2, the general solution is a linear combination (no quadrature is needed) 
x{t) = kiXi{t) + k2X2{t). A less known fact is the existence of other systems 
for which one can write the general solution as a (maybe nonlinear) function of a 
fundamental set of solutions and a set of constants determining each particular solution. 
The characterization of such systems is due to Lie [2] and these Lie systems appear 
quite often not only in mathematics but also in both classical and quantum physics 
jH El El [71 El [9l [TOl HI], the almost ubiquitous Riccati equation, 

^ = b2{t) x\t) + h,{t) x{t) + h,{t) . (5) 

which appears in the reduction process from a linear second-order equation when taking 
into account its dilation invariance, or in factorization methods, being one important 
example [11 El [121 US]- For such Riccati equation, if a particular solution is known, the 
general solution is obtained by two quadratures; if two particular solutions are known 
the problem reduces to one quadrature; and finally, when three particular solutions are 
known, Xi,X2 and X3, the general solution can be found from the cross ratio relation 

m 

X — Xi X3 — Xi 
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X - X2 X3- X2 
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which provides us a nonhnear superposition rule. 

All these properties can be better understood in the framework of Lie systems 
[2], conveniently extended in some cases to include second-order differential equations 
systems. These systems have a lot of applications not only in mathematics but also in 
many different branches of physics, both at the classical and the quantum level [3], some 
of them determined by t-dependent Hamiltonian that in general are not easy to handle. 
Furthermore, Lie systems are also relevant in Control Theory. 

The paper is organized as follows. Next section is devoted to first recall some 
important properties of Lie systems, and then to extend the theory for systems of second- 
order differential equations. Moreover, we illustrate the theory with some important 
examples with applications in physics and other relevant examples as Milne-Pinney 
[14] equation x = —uj'^{t)x + k/x^, which is usually studied together with the time- 
dependent harmonic oscillator y + oj^{t)y = 0, the so-called Ermakov system. Pinney 
showed in a short paper [15] that the general solution of the first equation can be written 
as a nonlinear superposition of two solutions of the associated harmonic oscillator. A 
generalisation of Ermakov systems and the corresponding Ermakov invariants, which 
appear as first integrals in a natural way, will also be analysed. It will be shown 
that all these properties can be better understood in the framework of Lie systems 
[2], conveniently extended in some cases to include second-order differential equations 
systems. The perspective of such systems as Lie systems is very important because 
they have a lot of applications not only in mathematics but also in many different 
branches of physics, both at the classical and the quantum level [3] . In particular, some 
t-dependent Hamiltonian systems are of this class, and Lie systems are also relevant in 
Control Theory. 

Section 3 is devoted to explain the reduction method which is the main ingredient 
in the setting of sufficient integrability conditions that are proposed in Section 4. This 
is based on the affine action of the group of curves in the associated Lie group on the set 
of Lie systems with such a group (see [16] for a geometric interpretation) and generalizes 
the results obtained in [1] for Riccati equation. The integrability conditions are used in 
Section 5 for studying Caldirola-Kanai oscillator [T71 118] and the possibihty of reduction 
to an autonomous system is proved as a consequence of the theory. Arbitrary time- 
dependent systems are also studied and explicit time-dependences are proved to satisfy 
the above mentioned integrability conditions and therefore to correspond to explicitly 
solvable models. 

2. Lie systems: A geometric approach 

The conditions for the system ([3]) to admit a superposition rule, i.e. there exists a 
function $ : M"(™-+i) — >■ M", x = . . . , ^i, • • • , such that its general 

solution can be written as x(t) = . . . , X(^rn)it); ki, . . . , kn), where | a = 

1, . . . , m} is a fundamental set of particular solutions and k = [ki, . . . , k^) is a set of 
n arbitrary constants, were studied by Lie [2]. Under very general conditions [5] the 
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necessary and sufficient conditions is that the associated t-dependent vector field X{x, t) 
can be written ClS db hnear combination 

r 

X{x,t) = Y,h^{t)X^^){x) , (6) 

a=l 

such that the vector fields {X(^a) \ a = 1, . . . ,r} generate a finite-dimensional real Lie 
algebra, i.e. there exist real numbers Cq,/?''', such that 

r 

[^(«)'^(/3)] = X^c«/3^^(7) ' a,/3 = l,...,r. (7) 

7=1 

Linear systems are particular examples whose associated Lie algebra is 0[(?t,,M) and 
m = n in the homogeneous case, or the corresponding affine algebra and m = n + 1 
in the inhomogeneous one. Riccati equation is another example for which X{x,t) = 
bo{t)Yo{x) + bi{t)Yi{x) + b2{t)Y2{x) with 

>^o(x) = ^, Y,{x)=x-^, Y2{x)=x'-^. (8) 
with commutation relations 

[Yo, Fi] = Yo, [Y2, Yi] = -Y2, [Y^, Yo] = -2Y,, 

closing on a s[(2, M) Lie algebra P US]. There is an action $Rie : SL{2, M) x M ^ M of 
the Lie group SL{2, M) on R = M U cxd given by 

^ r A \ C(X + (3 ,5 
$Ric(A,x) = — , if xjt — , 

"yx + 7 

$Ric(v4, oo) = a/7 , $Ric(A, -5/7) = 00 , , 
where the matrix A G SL{2,M.) is 

A = 

such that the fundamental vector fields of this action are the linear combinations of the 
vector fields Yq, Yi and 1^2- 

Another very relevant example of Lie system is given by a t-dependent right- 
invariant vector field in a Lie group G. A right-invariant vector field X^ is one such 
that X^{g) = Rg^:eX^{e). If {ai, . . . , a^} is a basis of TgG and define the right-invariant 
vector fields X^ by X^{g) = Rg^e^^a, then the t-dependent right-invariant vector field 
X{g,t) = —J2a=i^a(t)X^{g) defines a Lie system in G whose integral curves are 
solutions of the system g = — X]a=i ^ait) X^{g), and when applying Rg-i to both sides 
we see that g{t) satisfies 

r 

Rg-^{t)*g{t)git) = - ^ba{t)&a G T^G . (9) 
a=l 

Right-invariance means that it is enough to know one solution, for instance the one 
starting from the neutral element g{t), to know all the solutions of the equation with 
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any initial condition, i.e. we obtain the solution g'{t) with initial condition g'{0) = Qq as 
i?p(o)fi'(t). A generalization of the method by Wei- Norman [19] is very useful in solving 
such equation, and furthermore reduction techniques can also be used [20]. Finally, 
as right-invariant vector fields project on the fundamental vector fields in each 
homogeneous space for G, the solution of ([9]) allows us to find the general solution for 
the corresponding Lie system in each homogeneous space. Conversely, the knowledge of 
particular solutions of the associated system in a homogeneous space gives us a method 
for reduction the problem to the corresponding isotopy group [20j . This equation is also 
important because any Lie system described by a t-dependent vector field in a manifold 
M like ([6]) where the vector fields are complete and close the same commutation relations 
as the basis {ai, . . . , a^} determines an action $LieG : G x M ^ M such that the vector 
field Xct is the fundamental vector field corresponding to a^, and moreover, the integral 
curves of the t-dependent vector field are obtained from the solutions of thr equation ([9]). 
More explicitly, the general solution of the given Lie system are x(t) = $LieG(5'(^)) ^^o), 
where xq is an initial condition of the solution and g{t) is the solution for iQ with 
g{0) = e. 

The search for the number m of solutions and the superposition function $ have 
recently been studied from a geometric perspective [5]. Essentially we should consider 
'diagonal prolongations' to M"('"+^), X(a;(o), . . . ,X(m),t), of the t-dependent vector field 
X{x, t) = X^(t, x) d/dx\ given by 

m 

X(X(o), . . . = ^Xa{x^a),t) , t G M, 

a=0 

where Xa{x(^a),t) = X\x(^a)it) d / dx^f^^y such that the extended system admits n 
independent integrals of motion, which define in an implicit way the superposition 
function. 

As another example of Lie system of physical relevance, consider the t-dependent 
harmonic oscillator described by a Hamiltonian 

H{x, p,t) = 1 + i m{t)uj^{t)x'^ , 
/ m[t) 2 

whose Hamilton equations are given by 

OH p 



X 



dp m{t) 

p = ——— = —m(t)uj (t) X . 
ox 

Its solutions are the integral curves of the t-dependent vector field 

X(x,p,t) = -l-p^-m{t)uj\t)x^. (10) 
m(t) ox op 

If we consider the set of vector fields 
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which close on a sl(2, M) Lie algebra with commutation relations 

[Xo,Xi] = Xq, [X2,Xi] = —X2, [X2,Xq\ = — 2X1, 

then the t-dependent vector field (fTOj) can be written as a linear combination with t- 
dependent coefficients X{-,t) = J2l,=Qba(t)Xa{-), with 

fcoW = ^, h{t) = 0, b2{t) = m{t)uj\t) . 

Remark that the vector fields of this example and those of Riccati equation close the 
same commutation relations. 

As an instance of t-dependent harmonic oscillator we next develop the one with 
damping term, with equation of motion 

— (mo x) + mo fxx + k X = , k = mocu^ , 
at 

admits a Hamiltonian description with a t-dependent Hamiltonian fTfi [T8] 

H = exp(-/it) + l-niQ exp(/it)t^ V = ^ + ^m(t)t^ V (11) 
IniQ I 2m[t) 2 

with m(t) = mo exp(yut). The time-dependent coefficients as a Lie system are: 

^ exp{-fit) ^ _ , b2{t) = mo expifit) 

mo 

We can also consider the time-dependent frequency harmonic oscillator 

H=Ip' + lFit)ujy , (12) 

whose associated t-dependent vector field X is a linear combination X(-,t) = 
F{t)ujQX2{-) + Xo{-) with coefficients: 

6o(t) = l, 6i(t) = 0, b2it) = F{t)ujl 

The theory of Lie systems can be extended to include systems of second-order 
differential equations. Such a system whose associated system of first-order differential 
equations is a Lie system is called SODE Lie system [21]. This is the case of the 1- 
dimensional harmonic oscillator ([2]), where X{x,v,t) is a linear combination X(-,t) = 
X2(-) — ci;^(t)Xi(-) of the vector fields Xi{x,v) = xd/dv and X2{x,v) = vd/dx that 
close a finite-dimensional Lie algebra with Xs{x,v) = l/2{x d/dx — vd/dv). Similarly 
for the isotropic 2-dimensional case, 

= -u'^{t)xi 

= -Uj\t)X2 

but with Xi(xi, X2, fi, V2) = xi d/dvi+X2 d/dv2, X2(xi, X2, vi, V2) = f 1 d/dxi + V2 d/dx2 
and X3{xi,X2, fi, ^2) = l/2(a:i d/dxi + X2 d/dx2 — fi d/dvi — V2 d/dv2), in both cases 
with associated Lie algebra s[(2,M), 
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The general theory of Lie systems can now be used and so we can obtain a first 
integral, 0:2, fi, ^2) = X1V2 — X2V1, which is nothing but the Wronskian of two 

solutions, or the angular momentum from the physical viewpoint. This may be seen 
as a partial superposition rule: if Xi(t) is a solution of the first equation, then we 
obtain for each real number k the first-order differential equation for the variable X2, 
Xi{t) dx2/dt = k + Xi{t)x2, from where X2 can be found to be given by 

/■* dC 

X2{t) = k'xi{t) + kxi{t) ' ^ 



In a similar way, with three copies of the same harmonic oscillator, one obtains two first 
integrals, Fi{xi,X2,x,vi,V2,v) = xvi — xiv and F2{xi,X2,x,vi,V2,v) = XV2 — X2V that 
produce a superposition rule, because from 

xf 2 — X2V = ki 

XiV — ViX = k2 

we obtain the expected superposition rule [2T] : 

X = kiXi + k2X2 , V = kiVi + k2V2 . 

It has only recently been pointed out [21j that the differential equation [m [15] 



k 



X = -uj\t)x + , (14) 

X'^ 

is also a Lie system with associated Lie algebra s[(2, M): the corresponding t-dependent 
vector field is 

X{x, V, t) = v-^+ (-iu\t)x + ^]-^. (15) 



dx \ x^ J dv 

which is a linear combination X = L2 — uj'^{t)Li with 

d ^ k d d 

ov x'^ ov ox 

which are such that 

[Li,L2] = 2L3, [L3,Li]=Li, [L^,L2]-- 

with 

1 f d d 
L3 = - X— V — 

2 \ ox Ov 

Another important example of Lie system with the same Lie algebra s[(2, M) is the 
Ermakov system: 

X = Vx 
Vx = —Uj'^{t)x 

y = vy (16) 

1 



Vy = -uj^{t)y+^^ 



yC 
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with associated t-dependent vector field 

which is a hnear combination with t-dependent coefficients, X = —uj'^(t)Xi + X2, with 

d d d did 



dvx dvy dx dy y^ dvy 

This system is made up by two Lie systems closing on a sl{2, M) algebra: the ffist one 
is a harmonic oscillator and the second one is a Pinney equation. The system admits a 
ffist integral of the motion which can be found as a solution of XiF = X2F = because 
these two conditions imply that X^F = l/2[Xi,X2]F = 0. But XiF = means that 
F{x, y, Vx, Vy) = F{x, y,^) with ^ = xVy — yv^, and then X2F = is written 

^ OF ^ dF X dF _^ 
^ dx ^ dy y^ d^ 

and from the associated characteristics system we obtain 

xdy-ydx ^ y^ d^ d{x/y) _^ yd^ ^ ^ 
^ X ^ x 

the following first integral is found: 

ij{x, y, v^, Vy) = (^^^ = (^^^ + {xVy - yv^f 

which is the well-known Ermakov invariant. A possible generalisation of the Ermakov 
system ( IT6|l is given by: 

X = \f{y/x)-uj'^{t)x 

f , (17) 

y = -^9{y/x) -uj^{t)y 

y.i 

which reduces to f lTB]) for f{u) = and g{u) = 1. 

This system can be written as a ffist-order one by doubling the number of degrees 
of freedom by introducing the new variables Vx and Vy-. 

X = Vx 

Vx = -uj^{t)x + j^f{y/x) 

y = Vy 

. = -^^{t)y + jsgiy/x) 

which determines the integral curves of the vector field 



Integrability of Lie systems and some of its applications in physics 



9 



Such vector field can be written as a finear combination 

X = N2-(j\t)Ni 
where A^i and N2 are the vector fields 



9 d 
Ni = X— + y 



dv^r 



dVy 



d 1 d d 1 d 



dVy 



Note that these vector fields generate a three-dimensional real Lie algebra with a 
third generator 

1/(9 d d d 



2 Vdx 



dvx oy " dvy 



In fact, as 

[TVi, 7V2] = 27V3, [7V3, N,] = TV,, [Ns, N2] = -N2 

they generate a Lie algebra isomorphic to sl(2, M). Therefore the system is a Lie system. 

There exists a first integral for the motion, F : — > R, for any which can 

be found in a similar way and we will arrive to the first integral: 



1, ,2 r^y 

-{xvy-yv^} + j 



u 



ug \ - 

u 



du . 



This first integral allows us to determine a solution of one subsystem in terms of a 
solution of the other equation. 

Finally, we can revisit Pinney equation by considering a system made up by a 
Pinney equation in the x variable with two associated harmonic oscillator equations for 
variables y and z. by considering the system of first-order differential equations: 



X 


= Vx 


y 


= Vy 


z 


= 


Vx 


— —U) 


Vy 


— —U) 


Vz 


— —UJ 



k 

-(jj^{t)x+ — 

X'^ 



which corresponds to the vector field 



d_ 

' dx 



d 



d 



X = Vx— + Vy— Vz— + 



k d 



dz dv^ 



.,, , , (9 d d 



dv^ 



dv. 



dv. 



The vector field X can be expressed as X = A^2 — u;'^{t)Ni where the vector fields Ni 
and N2 are: 
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These vector fields generate a 3- dimensional real Lie algebra with A^3 given by 

^ 2 \ dx dvx ^ dy dvy dz dvz 
In fact, as 

[Nu N,] = 2iV3, [iVs, iVi] = iVi, [iV3, N2] = -N2 

they generate a Lie algebra isomorphic to sl(2, M). Thus this is a Lie system. 

There exist three first-integrals for the distribution generated by these fundamental 
vector fields: the Ermakov invariant Ii of the subsystem involving variables x and y, 
the Ermakov invariant I2 of the subsystem involving variables x and z, and finally, 
the Wronskian W of the subsystem involving variables y and z. They are given by 
W = yvz- zvy, 

= ^ (^{yvx - xvy)'^ + k ^ , ^2 = ^ (^{xv, - zv^f + k 

In terms of these three integrals we can obtain an explicit expression of x in terms 
of y, z and the integrals /i, I2, W: 

,/2 / , \ 1/2 

^ = V (^'^' ^ ^ \/4A/2 - kW^ yz 



We can recover in this way the result of [IS]: the general solution of Pinney equation 
can be expressed in terms of two solution of the corresponding harmonic oscillator 
problem with time- dependent frequency 



3. The reduction method 



Given an equation ([9]) on a Lie group, it may happen that the only non-vanishing 
coefficients are those corresponding to a subalgebra f) of g and then the equation reduces 
to a simpler equation on a subgroup, involving less coordinates. An important result is 
that if we know a particular solution of the problem associated in a homogeneous space, 
the original solution reduces to one on the isotopy subgroup pO]. 

One can show that there is an action of the group Q of curves in G on the set 
of right-invariant Lie systems in G (see e.g. [16] for a geometric justification), and we 
can take advantage of such an action for transforming a given Lie system into another 
simpler one. 

So, if g{t) is a solution of the given Lie system satisfying ((91) and we choose a curve 
g\t) in the group G, and define a curve g{t) by 'g{t) = g'(t)g{t), then the new curve in 
G, g{t), determines a new Lie system. Indeed, 

r 

Rm-^*9{t)(9it)) = Rg'~Ht)*9'{t)(9{t)) - ^6„(t)Ad(^'(t))a„ , 

0=1 
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which is similar to the original one, with a different right-hand side. Therefore, the 
aim is to choose the curve g'{t) in such a way that the new equation be simpler. For 
instance, we can choose a subgroup H and look for a choice of g'(t) such that the right 
hand side lies in T^H, and hence if ^(0) = e then 'g{t) G H for all t. This can be done 
when we know a solution of the associated Lie system in G/H, what allows us to reduce 
the problem to one in the subgroup H, see [20] . 

Theorem 1 Each solution of Ij^ on the group G can he written in the form g{t) = 
gi(t) h(t), where gi(t) is a curve on G projecting onto a solution gi(t) for the left action 
X of G on the homogeneous space G/H and h{t) is a solution of an equation hut for the 
suhgroup H , given explicitly hy 

{Rn-Uhh){t) = -AA{g^\t)) ^X^6a(t)a« + (i?^-i,^^(7i)(t)j eT,H . 

This fact is very important because one can show that Lie systems associated with 
solvable Lie algebras are solvable by quadratures and therefore, given a Lie system with 
an arbitrary G having a solvable subgroup, we should look for a possible transformation 
from the original system to one which reduces to the subalgebra and therefore integrable 
by quadratures. 

The result of the preceding Theorem proves that for a solvable Lie subgroup H of G 
there always exists a curve in G that transforms the initial Lie system into a new one in 
subgroup H of G. Nevertheless, it can be difficult to find out a solution of the equation 
in G/H that determines this transformation. Then, to be able to obtain one is more 
interesting to suppose also that this transformation is a curve in a certain subset of G, 
i.e. a one-dimensional Lie subgroup. When such a transformation exists, it is easier to 
obtain it, but it may be that such a transformation does not exist, and the conditions 
for the existence of such a transformation of a certain form are integrabihty conditions 
for the system. 

We could choose for showing this assertion a particular example: Riccati equation. 
One can find in the literature a lot of integrability criteria for Riccati equation, all of 
them particular examples of the above method. 

We can also consider more relevant examples in Physics, for instance, time- 
dependent harmonic oscillators. 

The results obtained for one system are valid for the other; they are essentially 
conditions for the equation in the group, and both are examples of Lie systems associated 
with the same Lie group: 5^(2, M). 

4. Integrability criteria for Lie systems 

Consider the particular case of Lie systems with associated Lie group SL{2, M), therefore 
valid for Riccati equation (I5l), Milne-Pinney equation ( |T4l) and the time-dependent 
harmonic oscillator described by (fTOj) . The group 5*^(2, M) contains the afiine group 
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(either the one generated by Xq and Xi or the one generated by Xi and X2), which is 
solvable. Therefore, a transformation from the given equation to one of this subgroup 
allows us to express the general solution in terms of quadratures. This happens when 
we know a particular solution xi of the given equation: x = Xi + z, what corresponds 
to choose 

1 —Xi 

1 



reduces the equation to 



dz 

— = (2a2Xi + ai)z + 02 z^ 
at 



The reduction by the knowledge of two or three quadratures has also been studied from 
this perspective and similarly Strelchenya criterion l22j. 

Each Riccati equation can be considered as a curve in and we can transform 
every function in M, under an element of the group Q of smooth S'L(2, M)-valued 
curves Map(M, ^L(2,M)), as follows: 

-i{t)x[t) + 5{t) -lit) 
e(Aoo) = a{t)h{t), e(A-5W/7W) = oo, 

when A = \ ) i ,) i \ eG . 

The image x'(t) = Q{A(t),x{t)) of a curve x{t) solution of the given Riccati equation 
satisfies a new Riccati equation with coefficients b[, 6q, given by 

62 = Ph2- 57 hi + 7^ 60 + 7*^ - (^7 , 

b'-^ = -2 (36 62 + {a6 + /?7) 61 - 2 07 60 + 5a - + /?7 - 7/3 , 
60 = 13"^ 1)2 — ajS bi + a"^ bf) + ajS — /3a . 

This expression defines an affine action of the group Q on the set of Riccati equations. 
As indicated above the vector fields defining the TDHO and the Riccati equation close 
the same conmutation relations and both examples are related with exactly the same 
equation in SL{2,M.), but they correspond to different actions. 

2 

Lie systems in SL{2, M) defined by a constant curve, a{t) = ^^CQ,aa, are integrable 
and the same happens for curves of the form a(t) = D{t) i ^c«a„ , where D is an 

\a=0 / 

arbitrary function, because a time parametrisation reduces the problem to the previous 
one. 

The system is essentially a Lie system on a one-dimensional Lie group. 

A straightforward application of the reduction method leads to the following result 
which is valid not only for Riccati equation but also for any other Lie system with Lie 
group S'L(2,M) in which the t-dependent vector field be described in terms of a set of 
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vector fields closing on the same commutation relations as those of Riccati equation. 
In this way, all these systems with the same t-dependent coefficients are related with 
the same equation in SL{2,M.) and this allows us to generalise for all these systems 
the integrability condition obtained in the next theorem. The other results can be 
generalised by means of changing the action of the Riccati equation for the action of 
the new Lie system. 

Theorem 2 The necessary and sufficient conditions for the existence of a transforma- 
tion: 

y' = G{t)y with G{t) > 
generated by the transformation 



= $Ric(A(t),|/) with A{t) 



and relating the Riccati equation in an interval I e 



a{t) 
a-\t) 



^ = bo{t) + b,{t)y + b2{t)y' , ytel, {bob2m ^ 0, 



with an integrable one given by 

D{t){co + Ciy' + C2y'^), 



dt 

where ci are real numbers, Ci e IR, is that for such constants cqC2 ^ and 

D {t)coC2 = bo{t)b2{t), = ci. 

The transformation is then uniquely defined by: 



b2{t)co 

As a consequence of this theorem, given the above Riccati equation if there are constants 
K and L such that 




boitMt) 

then there exists a time-dependent linear change of variables transforming the given 
equation into the solvable Riccati equation: 

'^^ = D{f:){co + cW + C2y'') 

where Ci = K, CoC2 = and D{t) is given as above. The existence of such constants 
K and L can be considered a sufficient condition for integrability of the given Riccati 
equation. Indeed, L is just used to make the square to exists and can always be obtained, 
the real integrability condition is the existence of the constant K. 
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5. Some applications in physics 

5.1. The Caldirola-Kanai oscillator 

Coming back to the time-dependent harmonic oscillator, we see that the Caldirola-Kanai 
model is an example of application of the previous Theorem. 

We have seen that it is a Lie system with group SL{2, M) and associated coefficients 

bolt) = — exp(-//i) , bi{t) = , b2{t) = moOUQ exp(//i) . 
mo 

Therefore, as bo{t)b2{t) = cUq, bi = and (62/62—^0/^0 = 2/x, wc sec that the integrability 
condition of the preceding theorem holds with K — /ji/ouo and L — 1 , and the function 
D is then a constant D = uq. Therefore we choose cq = l,ci = ///o^o and C2 = 1 and 
this example reduces to the system 





which can be easily integrated. This shows once again that the Caldirola-Kanai model 
can be reduced to an autonomous system, as a particular case of a more general situation. 

5.2. The t- dependent frequency harmonic oscillator 

In the case of the i-dependent frequency harmonic oscillator for which u;'^{t) — F{t)uJo, 
X{-,t) = F{t)uiX2{-)+Xo{-), i.e. 6o(0 = 1, hit) = and b2it) = F{t)ul the condition 
on F to satisfy the compatibility condition is 

I^^KujoVf 
and therefore F must be of the form 

F{t) 



{-Kuot + K'f ' 

for another constant K' and then the Hamiltonian which can be exactly integrated is 

2 2 {-KuJot + K'Y 
and the corresponding Hamilton equations are 



X — p 



^ i-Kuot + K')^'' 
Had we used instead the family of curves in SL{2,W) given by 

Ao{t) = 
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where Ui is a constant, and we would obtain the following relations among coefficients: 



V%2 + uiVbi + ulho - UiV 



b[ =bi + 2 



Ui 



V 



V V 



b' 



Assume we want to relate the t-dependent vector field X{-,t) = Xq + F{t)ujQX2, 
characterised by bo = 1, bi = 0, 62 = F(t)uQ, with another one characterised by 
6q, b[ and which is integrable, or more explicitly, related with the t-dependent vector 
field X{-, t) = D{t){coXo + C2X2), i.e. b'^ = Dcq, b[ = 0, and b'2 = Dc2, 

Then bi — b[ — 0, and the second equation shows that V — Ui, i.e. V{t) — Uit + Uo 
with uo e M, and using this condition in the first equation, together with bo — 1, 
it becomes 62 = V'^b2, and then, as the third equation gives us the value of D as 
CqD = Uq = l/V^, we see that b'^ — Dc2 — V'^F{t)ujQ, and therefore F must be 
proportional to {uit + Mo)~^: 



k 



C-2 



CoLol 



We can fix C2 = cOq and Cq = 1. Thus F{t) — {uit + Uo)~^. 

Finally, note that the time-dependent transformation Ao{t) which performs the 
reduction is 

X' 



X 



p 



V{t) 
— —uix + V{t)p 



which transforms the initial system of differential equations into 

1 . 




1 



p 



-u^x') 



Now, using the time-reparametrization 

r{t) 



dt' 



V^it') 



we obtain the following time-independent linear system 

dx' 



dr 
dp 



—LUqX 



whose general solution is 




sin(a;()T) 



cos(a;o''") 
—LJo sin(a;o''") cos(a;o'?") 
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and thus we obtain that 



x{t) = V{t) (cos(a;or(t))-^ + — sm{uJoT{t)){-uiXo + V{0)po) 
\ V [0) uo 

where we have used 

x'{0) = , p'{0) = -u,xo + V{0)po . 

The same computations are vahd for the Pinncy equation when F{t)Ljl is replaced 
by a time-dependent angular frequency u!^{t). 

Acknowledgments 

Partial financial support by research projects MTM2006-10531 and E24/1 (DGA) and 
a F.P.U. grant from Ministerio de Educacion y Ciencia are acknowledged. 

References 

[1] J.F. Carificna and A. Ramos, Int. J. Mod. Phys. A 14, (1999) 1935- 51. 

[2] S. Lie, Vorlesungen iiber continuierliche Gruppen mit Geometrischen und anderen Anwendungen, 

Edited and revised by G. Scheffers, Teubner, Leipzig, 1893. 
[3] P. Wintcrnitz, Lie groups and solutions of nonlinear differential equations, in: Nonlinear 

Phenomena, K.B. Wolf Ed., Lecture Notes in Physics 189, Springer- Verlag, N.Y., 1983 
[4] J.F. Carinena, J. Grabowski and G. Marmo, Lie-SchefFers systems: a geometric approach, 

Bibliopolis, Napoli, 2000. 
[5] J.F. Carifiena, J. Grabowski and G. Marmo, Rep. Math. Phys. 60, (2007) 237-58 
[6] R.L. Anderson, J. Harnad and P. Winternitz, Lett, in Math. Phys. 5, (1981) 143-148. 
[7] J. Harnad, P. Winternitz and R.L. Anderson, J. Math. Phys. 24, (1983) 1062-72. 
[8] S. Shnidcr and P. Winternitz, J. Math. Phys. 25, (1984) 3155-65. 
[9] M.A. del Olmo, M.A. Rodriguez and P. Winternitz, J. Math. Phys. 27, (1986) 14-23 
[10] M.A. del Olmo, M.A. Rodriguez and P. Winternitz, Integrability, chaos and nonlinear superposition 
formulas for differential matrix Riccati equations, in: Quantum Chaos and Statistical Nuclear 
Physics, pp. 372-378, Lecture Notes in Physics 263, 1986, 
[11] M.A. del Olmo, M.A. Rodriguez and P. Winternitz, J. Math. Phys. 28, (1987) 530-535. 
[12] J.F. Carinena and A. Ramos, Acta Appl. Math. 70, (2002) 43-69. 
[13] J.F. Carifiena, A. Ramos and J. do Lucas, Elect. J. Diff. Eqs. 122, (2007) 1-14. 
[14] W.E. Milne, Phys. Rev. 35, (1930) 863-67. 
[15] E. Pinney, Proc. A.M.S. 1, (1950) 681. 

[16] J. F. Carinena and A. Ramos, Lie systems and Connections in fibre bundles: Applications in 
Quantum Mechanics, in: S**" Int. Conf. Diff. Geom and Appl, p. 437-452 (2004), J. Bures et al 
Eds., (Matfyzpress, Praga, 2005). 

[17] P. Caldirola, Nuovo Cim. 18, (1941) 393-400. 

[18] E. Kanai, Prog. Thcor. Phys. 3, (1948) 440 42. 

[19] J. Wei and E. Norman, J. Math. Phys. 4, (1963) 575 81. 

[20] J.F. Carifiena, J. Grabowski and A. Ramos, Acta Appl. Math. 66, (2001) 67-87. 

[21] J.F. Carifiena, J. de Lucas and M.F. Rafiada, Nonlinear superpositions and Ermakov systems, in: 

Differential Geometric Methods in Mechanics and Field Theory, pp. 15-33, eds F. Cantrijn, M. 

Crampin and B. Langerock (Academia Press, 2007) 
[22] V. M. Strelchenya, J. Phys. A: Math. Gen. 24, (1991) 4965-4967. 



